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Abstract We investigate the effects of low energy heavy nuclei 
from the supernova explosion on nearby interstellar space. In addition 
to the ionization ahd heating of the Gum Nebula, these particles may 
produce detectable fluxes of X-rays and gamma rays, both as continuum 
radiation and line emission. 


Paper presented by R, Ramaty at the Symposium on The Gum Nebula and 
Related Problems, Goddard Space Flight Center, Greenbelt, Maryland, 

May 19, 1971. To appear in The Gum Nebula and Related Problems , edited 
by S. P. Maran, J. C. Brandt and T, P. Stecher. 
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INTROUUCTTON 

The possibility that; the ionization of the Gum Nebula may have 
been produced by energetic charged particles from the supernova Vela X has 
been suggested by Ramaty et al. (1971). The principal features of this 
model are the existence of a prolonged source of ionization and heating 
and the possibility of direct detection of hard photon emissions from 
contemporary fluxes of energetic particles in the nebula. In the present 
paper, we shall summarize the model and discuss its main features. We 
shall also present an estimate of gamma-ray line emission from the inter- 
action of energetic nuclei with the ambient gas and we shall discuss 
some additional ideas regarding the maintenance of a possible high 
temperature in the filaments of the nebula. 

Particle Propagation and Energy Deposition 

A necessary requirement of a cosmic-ray model of the Gum Nebula is 
that in a time less than or equal to the age of the supernova remnant 
the particles should be capable of both reaching the outer edges of the 
nebula and depositing a significant fraction of their energy as ionization 
loss to the ambient gas. 

A lower limit on the velocity eg of the ionizing nuclei can be 
obtained by assuming rectilinear propagation of particles in the nebula, i.e. 

p > L/ct (1) 

where L and t are the nebular radius and age respectively. If we adopt 
an age of about 10^ years, as deduced from the period and rate of change 
of period of PSR 0833 (Reichley, Downs and Morris, 1970), and a nebular 
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radius of about 360 pc, as obtained (Alexander et al,, 1971) from Lya 

absorption data and a local hydrogen density of 0,24 cm”^, the lower limit 

on 8 is about 0.1 corresponding to an energy per nucleon 6^5 MeV/nucleon, 

It should be remembered that this lower limit is sub ject to considerable 

uncertainties in both t and L, In particular, if the age is larger than 
4 

10 years (Shklovsky, 1970), the Gum Nebula may have been produced by 
even lower energy particles. 

Cosmic-ray particles in general do not propagate along rectilinear 
paths but follow magnetic field lines and are scattered by irregularities 
in the field. The mean distance traversed in such diffusive motion may be 
approximated by 

where X is the diffusion mean free path. For L ^ 360 pc and t 10^ years, 

XB =“ 20 pc. If 8 = 0.15 (£ = 10 MeV/nucleon), X =“ 130 pc. This value 
is quite large but not inconsistent with recent ideas on cosmic ray 
propagation (Ramaty and Lingenfelter , 1971). In particular, if cosmic rays 
propagate by compound diffusion (Lingenfelter, Ramaty and Fisk, 1971) or 
if cosmic rays are confined to the galactic disk by reflecting boundaries, 
both the low anisotropy and abundances of fragmentation products of 
galactic cosmic rays are not Inconsistent with the long diffusion mean 
free path required for the ionization of the Ginn Nebula. Finally, a 
relatively scatter-free motion in the nebula is also supported by the low 
values of both the ambient gas density and magnetic field in this region 
of interstellar space. 
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The ionization loss rate and the time for the deposition of halt 
the energy of a nonrelativistic particle of mass number A and nuclear 
charge Ze in neutral hydrogen of density nn ^^e given by 

c: lo" MeV/nucf. s*c 

and ^ . 

A/('f.jjn.«) aears 

where £ and nj^ are in MeV/nucleon and cm"^, respectively. As a result 
of electron capture during the slowing down process the effective charge 
Zeff becomes smaller than the nuclear charge Z, Zeff/Z may be approximated 
by (Pierce and Blann, 1968) 


At 




Since even for Fe^^ at 10 MeV/nucleon Zeff/Z=0.9, we shall use Z instead 
of Zeff in the evaluation of equations (3) and (4), 

In Table 1, given for various 10 MeV/nucleon particles. 

As can be seen, in 10^ years these particles will deposit a major fraction 
of their energy only if n^^ » 1 cm"^. The average electron density < ne> 
in the Gum Nebula is only 0.16 cm"^, but the density in filaments is much 
higher, ne X < ne > where X 65 in the clumping factor in the nebula 
(Brandt et al,, 1971), A possible model, therefore, is one in which initially 
the nebular region consisted of neutral clouds or clumps with hydrogen densities 
approximately equal to the presently observed electron densities in filaments. 
The particles ejected from the supernova propagate in a scatter- free manner 
in the intercloud medium and subsequently penetrate clouds and lose their 
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energy in a tine scale shorter than the trapping time in the cloud, A 
possible trapping mechanisms, based on wave-particle interactions, has 
been discussed by Ramaty et al, (1971). 

If the initial hydrogen density in clouds is equal to the present 

_3 

electron density of about 10 cm , it is very unlikely that low energy 
protons have ionized the Gum Nebula, since, as can be seen from Table 1, 
it would take about 2x10^ years for these particles to lose half their energy. 
On the other hand, if the ionization were produced by heavy nuclei such 
as Si and Fe , these particles would deposit sizable fractions of their 
energy in clouds on a time scale of 2000 to 3000 years leaving about 7000 
to 8000 years for scatter-free propagation in the intercloud medium. This 
is consistent with the estimates of the transit times to the edges of 
the nebula given above, 

A specific consequence of this model is that the degree of ionization 
of the intercloud medium is much smaller than that of the filaments. This 

4 

is not the case for models in which the ionization is produced by photons 
since these lead to the complete ionization of both the clouds and inter- 
cloud medium out to an appropriate Stromgren radius. 

Mass and Energy Requirements of Supernova Ejecta 

A fast particle of energy ^0.1 MeV/nucleon expends some 36 eV to 

produce an ion pair (Bethe and Ashkin, 1953), Therefore, if the nebula 
6 2 

contains Ne 2x10 electrons, as suggested by Brandt et al, (1971), the 

52 

particles ejected by the supernova must lose at least 1,2x10 ergs. For 
the present value of Ne used, this is a lower limit, since a fraction of 
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the energetic particles may slow down In an Ionized medium thereby 
converting their energy directly to heat, and the total number of electrons 
produced could be larger than Ne since some recombination In the filaments 
may have already taken place. 

The velocity distribution of the ejected material from a supernova 


explosion may be approximated by (Colgate and McKee, 1969) 


XTr 


'’■i ( >-F) 


•/3 


-'h 


0.1S z F IFo * 3 /^ 

Fo 


( 6 ) 


io.67^F , 

where F Is the external mass fraction. The corresponding number spectrum 
of nuclei of mass number A In the nonrelatlvlstlc region Is given by 

'(f.ny'"' 


m 


7 Ai^f i 6 >€„ 


(7) 


Upon Integration, the total kinetic energy W of the ejected material Is 

\J - c Mej c 

^ - y is; 

where mp Is the proton mass. 

Since W > 1. 2xl0^^ergs , the ejected mass In units of a solar mass 
and the characteristic ejection energy In MeV/nucleon must obey 


— ^ £© ^ & M«.V/hMcfton (9) 

Mo 

As discussed above, £ cannot be smaller than about 5 MeV/nucleon, 

O 

since otherwise the particles will not reach the edges of the nebula in the 
available time of about 10^ years. On the other hand, cannot be much 
larger than about 20 MeV/nucleon, since, as can be seen from Table 1 and 
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equation (4), the ejected nuclei, even If they are all Fe^^, will not 

have enough time to deposit their energy In clouds. For the range 

5 < Sq< 20 MeV/nucleon, the lower limits on the ejected mass range between 

0.25 Mq to 1 Mg. For the efficient Ionization of the nebula, the ejected 

12 

material must consist of nuclei heavier than at least C . 

We now proceed to examine the observational consequences of the model. 


Hard Photon Emissions by Energetic Particles In the Nebula 

X-ray emission by bremsstrahlung Is a necessary consequence of the 
Interaction of fast nuclei with Interstellar gas (Boldt and Serlemltsos, 

1969). The photon energy spectrum that would result from the particle 
distribution given In equation (7) Is essentially flat with an end point 
hvQ at (mg/mp)£^, where mg and mp are the electron and proton masses, 
respectively. 

Let M/(Amp) nuclei of mass nymber A and atomic mnnber Z be currently 
trapped In clouds of density ne. The average bremsstrahlung volume emlsslvlty 
In the Gum Nebula (Independent of the Ionization state of the gas) Is given by 

r H V ^ ^<-f> 

-25 2 

where Qq = 6.7x10 cm , o = 1/137, V Is the volume of the Gum Nebula, and 

< c0 > Is the particle velocity averaged over the -inergy spectrum of the 

particles. For < 0 > * 0,15, Z^/A * 8.7 (average of Sl^® and Fe^^ In a 

-3 63 3 

ratio of 2 to 1), ne = 10 cm and V * 1,35x10 cm (Brandt et al, , 1971) 


we get 


SaIo 


-28 ^ 

lift 




cwT^ sea-' 


( 11 ) 
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According to the geometry of the Gum Nebula described by Brandt et 
al, (1971), we expect it to appear as a disk source of X-rays, extended 
in galactic longitude (AX 90°) and relatively thin in galactic latitude 
(A^ =“ 7°). Collecting radiation from an Interval of length L penetrating 
the nebula, the intensity viewed by a detector of aperture (A^)q < A^ 
would be expressed by 
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rtuL 

Using a detector of (A^)q = 4°, viewing at X (260 + 40)°, which 

spans the Gum Nebula, Cooke, Griffiths and Pounds (1969) observed an 

enhanced X-ray emission associated with the galactic disk of about 
-9 7 - 1 

3x10 ergs ( 6 m sec rad)” above 1,4 keV, with an essentially flat energy 
spectrum extending to about 10 keV. The upper limits to this flux at 
energies above 12.5 keV, as set by Hudson, Peterson, and Schwartz (1971), 
provide evidence that the spectrum is characterized by a break in the 
vicinity of about 10 keV.^ Using the volume emissivity as obtained from 


equation (11), L = 400 pc and (A^q= 4°, we get dl/dX = 6x10 ^ ergs 

-1 -2 -1 

sec cm rad . Since a significant part of the apparently diffuse 

background observed from the direction of the Gum Nebula may be coming 

from undetected discrete sources (Cooke and Pounds, 1971), *'he flux of 
-9 -2-1 -1 

3x10 ergs cm sec rad must be considered as an upper limit. This then 
sets an upper limit M/M 0 < 0.5 on the mass of energetic charged particles 
currently trapped in clouds. 

In addition to bremsstrahlung X-rays, contempoary fluxes of energetic 


heavy nuclei in dense clouds may also produce detectable gannna-ray line 
emission in the few MeV range. As an example, we consider the 1.78 MeV 


line of Si . The cross section Oy for the reaction Si^® (PiP*) 
has a threshold at about 5 MeV, peaks at about 9 MeV and has a value 
of about 300 mb at 10 MeV (McGowan et al,, 1969, 1970), 

As above we assume that M/ (Amp) particles are currently trapped in 
clouds. The average volume emissivity of the Gum Nebula in line emission 
at a photon energy is then given by 

where np equal to ne Is the proton density in clouds, and the factor of 
2/3 is the assumed fraction of silicon nuclei in the superru>va ejecta. 
For the values of np“ ne, V and < 8 > given above 


Since present measurements of gamma rays at these energies have been made 
with omnidirectional detectors only (Vette et al., 1970), we shall make the 
simplifying assumption that the flux at earth is produced by a point source 
at the center of the Gum Nebula. The flux seen by an omnidirectional 
detector will then be given by 




For r “ 460 pc, « 1.78 MeV and M/Mg < 0,5, F (1.78 MeV) < 8x10" photons 

cm"^sec~^. This flux should be compared with the omnidirectional counting 
“1 9 1 

rate of 1.2x10 photons cm ^sec" in the 1 to 2 MeV range (Vette et al,, 
1970). Since the energy of the photon-emitting nuclei is about 10 MeV/nucleon 
or less, line broadening should be less than 157.. Therefore, detectable 
ganraa-ray line emission could be observed with detectors of moderate energy 


resolution 
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It Is of some interest to speculate about line emission at 4.43 

MeV from and 6.14 MeV from 0^^. The cross sections for the excitation 

of these levels were plotted as a function of energy by Lingenfelter and 

Ramaty (1967). At 10 MeV/nucleon Oy 200 mb for and “* 100 mb for 

0 , If and Mq are the present masses of energetic caihcn and oxygen 

nuclei in clouds the expected photon fluxes from the Gum Nebula are 

F(4.43 MeV) =“ 4x10“^ photons cm^sec ^ and F(6.14 MeV) =“ 1,5x10”^ 

Mo/M^ photons cm“'^sec" . When compared with the observed counting rate of 
-2 -2 -1 

about 5.6x10 photons cm sec in the 4 to 6 MeV band (Vette et al., 1970), 
we see that line emission from carbon and oxygen could be detected and, 
moreover, it could account for at least part of the excess radiation above 
the intergalactic background as calculated by Brecher and Morrison (1969). 

In addition to the gamma ray region, line emission should also be 
produced in the X-ray region as a result of charge exchange between energetic 
nuclei and H atoms resulting in the capture of electrons to excited states. 
Cascades down to the ground state occur, which produce the analogue of 
Lya and LyB emission for the case of capture of K-electrons. Silk and 
Cteignan (1969) have given a discussion of th'" nrocess, with application 
to low energy galactic cosmic rays. The contemporary flux of X-ray line 
emission, however, is difficult to determine, since unlike bremsstrahlung 
and nuclear collision which do not depend on the ionization state of the 
ambient matter, charge exchange takes place in neutral hydrogen only. 
Estimates of the time integrated fluxes in the 7 keV line of iron and 
2 keV line of silicon were given by Ramaty et al. (1971). 
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Cooling of the Nebula 

Alexander et al. (1971) suggested that the temperature of the 
nebula as a whole is quite high, probably in excess of 3xlO‘^K. While it 
is possible that the dense filaments could at present be cooler than 
lO^K, it is of some interest to examine the effects of a contemporary source 
of ionization and heating on the thermal state of the filaments. 

According to Cox and Tucker (1969), for a gas of cosmic abundances 
in ionization equilibrium, the radiative loss rate at SxlO^K is P 2x10”^^ 
ng^ergs cm“^sec"^. The corresponding cooling time is tc"^ 3/2 ngkT/P =“ 

1700/ng years. For ng = 10 cm"^, tg 170 years, much shorter than the 
age of the nebula. 

The rapid cooling of the filaments can be partially compensated by a 
contemporary source of heat such as a large flux of energetic particles 
presently trapped in clouds. The rate of energy loss of a particle of 
mass number A and nuclear charge Ze in a plasma of density ne is approximately 
^.jiven by equation (3) w?.th n^j replaced by 4ne (Hayakawa and Kitao, 1956). 

If M/(Ainp) nuclei are currently trapped in clouds, the contemporary heat 
input is given by 


a S.e>f 

For Z^/A = 8.7, ^ = 10 MeV/nucleon, -ne = 10 cm"^ and M/Mq < 0.5, Q < 0.6P. 


Thus, even though a contemporary flux of energetic particles, equal to the 
upper limit determined from the calculated bremsstrahlung emission and the 
observed disk component of diffuse x-rays could make a contribution to the 
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malntenance of the teinperature In the filaments, It is not quite 

sufficient to balance the radiative energy loss of the plasma. The 

characteristic time scale for contemporary heating, ty,, is given by 

2 

equation (4) with np replaced by 4ne. For Z /A = 8,7, £ = 10 MeV/nucleon 
and ne = 10 cm”^, tj^ =*■ 800 years. 

The radiation loss P, however, could be smaller than the value 
calculated by Cox and Tucker (1969) if the abundances of carbon and 
oxygen in the Gum Nebula are below universal abundances or if the ions 
to and 0^^ to are significantly depleted from their values 
corresponding to ionization equilibrium by the initial ionization process. 

If the maintenance of a high temperature in the filaments is achieved by 
the latter mechanism, the characteristic time of interest is the recombination 
time of 0"*^ to and to 

Up to temperatures of a few times lO^K dielectronic recombination 
for 0 and C is negligible. Using Seaton's (1959) hydrogenic approximation, 
the radiative recombination coefficient to all levels of an ion X is 
given by (e,g, Allen and Dupree, 1969) 

where X = 157890 (mj^) /T, For ne = 10 and T = 2x10%, the recombination 
time for both O"^ -» O'*’^ and -♦ is about 900 years. By combining 
this time scale with the contemporary heating time tjj» we see that the onset 
of rapid cooling could be delayed by as much as 1500 to 2000 years, A high 
temperature in the filaments could then be consistent with a delayed 
mechanism for the ionization and heating of the Gum Nebula. 
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lonization and Heating of the Interstellar Medium 

Several recent theoretical descriptions (Pikel'ner 1967; Balasubrah- 
manyan et al, 1968; Spitzer and Tomasko 1968; Spitzer and Scott 1969; 
Goldsmith, Habing, and Field 1969) of ionization and thermal equilibria 
for interstellar HI regions, based upon current models for the interstellar 


gas, indicate that the rate of ionization per hydrogen atom is 


10 
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(sec H atom) . Hjellming, Gordon, and Gordon (1969) find that observed 

pulsar dispersion measures may be best fitted with such a model for = 
-15 -1 

(2,5 + 0,5) X 10 (sec H,'atom) . For an equilibrium situation, the 
rate of electron-ion recombination is the most direct measure of ; 
recent observations of H8 hydrogen line emission from interstellar HI 
(Reynolds 1971) yield a direct measure of the recombination rate, and the 
corresponding ionization rate in the regions examined could be as high as 
10"^^ (sec H atom) 

Since essentially all the gas in the Gum Nebula has been ionized by 
the supernova explosion, the average rate of ionization in interstellar 
hydrogen by such events would be giver, by 


( 18 ) 

where < n^^ > is the mean interstellar hydrogen density and fgj^ is the 

average supernova frequency per unit volume in the galaxy. If supernova 

66 3 

occur at a rate of 1 per 100 years in a volume of 4 x 10 cm"^, fgj^ 

8 X 10 cm"^sec~^. For Ne =“ 2 x 10^^ electrons (Brandt et al, 1971), 

1 A "X "3 

we obtain ^ 1,6x10”^^ sec cm , Radio observations at 21 cm 

indicate a mean value for the neutral hydrogen density in the galactic plane 
of about 0,7 cm (Spitzer 1968); however, optical depth effects probably 
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raise < nn > to ~ 1 cm , A further increase by about a factor of 2 in 
the value of < n^^ > may be due to the presence of molecular hydrogen in 
dense clouds. Hence the mean galactic value of as indicated by the 

parameters deduced for the Vela X supernova is about 8x10 sec”^ (H atom) \ 
This ionization rate could be reduced if the frequency of supernova explosions 
of the type that produced the Gum Nebula is lower than the assumed rate of 
1 per 100 years or if the actual free electron content of the Gum Nebula 
is less that that estimated by Brandt et al, (1971). Furthermore, even if 
the average supernova output is similar to that of Vela X but occurs in a 
region with larger magnetic fields, so that the particles propagate over a 
shorter distance before they lose their energy, the ionized volume produced 
by the supernova would be smaller than that of the Gum Nebula, In this case, 
the total ionization per supernova would be reduced since a larger fraction 
of the total energy would go into heating. Also, those supernovae which 
occur at more than about 100 pc away from the galactic plane (e,g. the Crab) 
will produce energetic particles which may escape from the disk without 
producing significant ionization. 

In order for HII regions such as the Gum Nebula to merge into the 

HI of the interstellar medium, the recombination and cooling times of the 

clouds should be much less than the time between supernova explosions tgN 

-77 -3 -1 

in the volume of the nebula. For f„„ = 8 x 10 qm sec and V = 1.35 x 

SN 

lO^^cm^, tgjj 3 X 10^ years. 

0.7 

The recombination time of ffydrogen is approximately t^.^^ 130 T * /ne 

(Bates and Dalgarno 1962), For ne = 10 cm ^ and T = 5xl0^K , t^.^^ ^ 2.5x10^ 
years. Since the recombination rate increases with decreasing temperature. 
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both t(. and tj-g c much snialler than so that the ionized clouds 

in the Gum Nebula will rapidly merge into the HI of the interstellar medium, 

-3 

In the intercloud medium, for a density < ne > 0,16 cm , the 

temperature will not decrease appreciably below lO^K in 3x10^ years. The 
recombination time then becomes comparable to tgjj so that the degree of 
ionization is approximately constant as a function of time. The fact that 
the inferred degree of ionization of the intercloud medium probably does 
not exceed about 10% (Hjellming, Gordon and Gordon, 1969) is consistent 
with a cosmic-ray model for the ionization of the Gum Nebula, since, as 
discussed above, the energetic charged particles from the supernova will 
deposit most of their energy in clouds and not in the intercloud medium. 
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